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We study the effects of sudden change in the sound velocity on primordial curvature perturbation 
' spectrum in inflationary cosmology, assuming that the background evolution satisfies the slow-roll 

t-H ' condition throughout. It is found that the power spectrum acquires oscillating features which are 

f*"^ , determined by the ratio of the sound speed before and after the transition and the wavenumeber 

£Nj ■ which crosses the sound horizon at the transition, and their analytic expression is given. In some 

— «. ' values of those parameters, the oscillating primordial power spectrum can better fit the observed 

' , Cosmic Microwave Background temperature anisotropy power spectrum than the simple power-law 

power spectrum, although introduction of such a new degree of freedom is not justified in the context 
of Akaike's Information Criterion. 

as ■ 

PACS numbers: 98.80.-k, 98.90.Cq 

o 

U: 

J3 I. INTRODUCTION 

Standard inflationary cosmology [H-Q predicts nearly scale-invariant power spectrum of the primordial perturbation 
■£3 ', 0-0] ■ Such a power-law like perturbation spectrum A 2 oc fe" 4-1 , where ( is the comoving curvature perturbation 
^ ■ and n s is the scalar spectral index with n s ~ 1 has been prefcred also from a number of obserbations Q. If we 
1— —i[ look into the detailed structure of the Cosmic Microwave Background (CMB) temperature anisotropy, some hints 
of small deviations from the simplest form of the curvature perturbation power spectrum show up. Among those, 
anomalously low values of the quadrupole moment or several sharp glitches in the large scale WMAP data corresponing 
to I ~ 20 — 40 are famous and most intensively studied Furthermore, several groups have reported strong evidence 
of the deviation from a simple power-law type power spectrum by reconstructing the primordial power spectrum from 
(*T) . the WMAP data (Tol-[l9|. In particular, Ichiki ct al. [lo, HH claim that they found an oscillatory modulation localized 
t^J" ' around the comoving wavenumbcr k ~ 0.009[Mpc _1 ] (£ ~ 120) in the power spectrum at 99.995% confidence level. 
, Theoretical calculation for the power spectrum has been also sophisticated recently and some inflation models based 
■ on a realistic high-energy physics c an g enerate peculiar features on the curvature perturbation. Those include the 
so-called Trans-Planckian effect [22H24I ]. particle production due to the coupling between the inflaton and another 
7— I • scalar filed [25T - [27j , temporal violation of the slow-rolling of the inflaton field |28rl35| , and some other models [HI, [I?} . 
^ ' To go further into the accurate cosmology, it is crucial to study further the details of the primordial perturbations. 
• i-h , This leads us to the detailed structure of the inflaton Lagrangian or the true high-energy physics that has realized in 
1 our Universe. 

In this paper, we concentrate on the role of the sound velocity, which is defined as the propagating speed of the 
■ - - > linear perturbation in the next section. In some high-energy physics theories, there appear non-canonical kinetic terms 
in the Lagrangian and in that case the sound velocity deviates from unity j3q.[39l|. Furthermore, if those kinetic terms 
of the inflaton field couple to some time-dependent variables which can be seen, for example, in DBI inflation scenario 
j40l | , the values of the sound velocity can change during inflation (4lTl43j . Then we have to consider the possibility of 
some non-trivial dynamics of the perturbation. As a result, new degrees of freedom such as the sound velocity may 
be observed in the CMB temperature fluctuation and be tested in the future high-precision CMB observations such 
as PLANCK or CMBpol. 
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This paper is organized as follows. In the next section, we comment on the background evolution in our scenario and 
introduce the sound velocity. In sec IIII1 the basic variables and its evolution equation for the curvature perturbation 
are described. Then we consider the paricular types of the variation of the sound velocity. The first type is a step-like 
function, which is discussed in Sec. IIV1 and the second type is a top-hat type function, which we will study in Sec. 
IVl The last section is devoted to the summury and discussion. 



II. BACKGROUND ASSUMPTION AND SOUND VELOCITY 



In the standard single inflaton field with a canonical kinetic term, the sound velocity c s defined by the propagation 
speed of linear perturbation has the value equal to the speed of light, namely, c s = 1. This is easily checked by 
considering the action of the canonical inflaton field, 

S = I d A xV~9 (\r + X + V (</>)) , X = -\g^d^d v ^ (2.1) 



2 v 7 2 

where R denote the Einstcin-Hilbcrt action and we set the reduced Planck scale to unity (8irG = 1). Expanding the 
action around the background Friedman-Robertson- Walker (FRW) metric up to second-order, then we can see that 
the sound velocity which appears as the coefficient of the spatial derivative term is exactly equal to one. 
Generalizing the above inflaton action to an arbitary function of the kinetic term X , 

P(X,<f>)Y (2.2) 

the situation changes completely. In this action, expanding around FRW metric up to second-order, we find the sound 
velocity is given by 



2P XX X + P x 



(2.3) 



where the subscript X represents a derivative with respect to X . The specific second-order action for the approapriate 
perturbation variable, what we call the comoving curvature perturbation, will be presented in the next section. It 
is clear from (|2.3p that, in principle, the sound velocity takes various values corresponding to the functional form 
P(X, 4>) and the values of X and <fi. 

In this paper, in order to extract the effects of the change of the sound velocity alone, we study the cases it changes 
suddenly once or twice during inflation without affecting the background evolution. The sudden change of the sound 
velocity can be realized, for example, if the Lagrangian has terms like 

P(A^) D /(0)X + [1 - /(0)]X 2 /A 4 , /(0)= ( l_, )/d , (2.4) 

or 

P(X,cj ) )Db(cj ) )X + X 2 /A 4 , b^) ~ {<j> - faf + 0(0 -&))"), (n>3) (2.5) 

where (/) is a scaler field which we assume to be the inflaton and A is some cut-off scale. In the former case, we find 
c 2 = 1 when (j> <C <fio and c 2 = 1/3 when <f> ^ tfio. Such a transition of the sound velocity due to the motion of the field 
(j> is approximately described by a step function. In the latter case, c 2 = 1/3 only when ~ O an d at the other values 
of </>, c 2 = 1. This variation is well approxmated by a top-hat type function. In generic unified theories including 
string theory, there appears many scalar fields which have non-trivial kinetic terms and it is likely that those kinetic 
terms naturally couple to the scalar field itself in a complicated way. Here we consider simple tractable examples of 
posssiblc deviation from the canonical model as described above. 

If we take these models as they are, the background evolution, in particular the slow-roll parameters, may also be 
severely modulated unless some fine-tunning is applied. Since the effects of the sudden change of slow-roll parameters 

He . . 

« = * = lTc> (2 ' 6) 

where a dot represents a derivative with respect to the physical time t, have been already studied in the literatures 
[3lL [3^1 . we concentrate on the cases only the sound velocity changes suddenly as mentioned above. Thus, we impose 
the slow-roll conditions 



|e|<l, M<1, 



(2.7) 
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throughout this paper. Such a situation can be also realized in the curvaton scenario as discussed in Appendix A. 

In the usual calculation of linear perturbation in this type of generalized non-canonical kinetic term inflation models, 
new parameters which parametrize the time variation of the sound velocity such as 

Vs = (2.8) 



c s H' " e s H 

are introduced and additional slow-roll conditions |e s | < l,|f? s | < 1 are imposed. In this paper, however, being 
interested in the effect of the more general variation of sound velocity mentioned above, we consider the situation 
where these conditions are temporarily violated by the sudden change of the sound velocity. 

In the following sections, we derive the approapriate variables for calculating the perturbation and the matching 
condition at the transition epoch. Then we evaluate the final power spectrum for two concrete examples of the sound 
velocity variation. 



III. CURVATURE PERTURBATION AND BASIC EQUATION 

The scalar perturbation from the Background FRW metric in the conformal Newtonian gauge is incorporated as 

ds 2 = a 2 [-(1 + 2^)dr 2 + (1 + 2$)5 ij dx i dx j ] (3.1) 

Denoting the inflaton perturbation 8<j) and the background value of the inflaton field </>, the most common perturbation 
variable, namely, comoving curvature perturbation ( is defined as j43 - l46j 



C = $ + ^50. (3.2) 

The basic action and the equation of motion in the linear theory for £ is written in terms of 

v = Cz, (3.3) 

where a new parameter z is defined by 

z^^. (3.4) 

C s 

The second-order action for v is derived by substituting (|3.ip and <fi = (f> + 5<fi into (|2.2[) , neglecting the higher-order 
terms and using the background evolution equation. The result is 

SW = \ [ drd 3 x ( v' 2 + c 2 s vAv + —v 2 ) , (3.5) 



2 J \ * 

where a prime denotes a derivative with respect to the conformal time r. Then, the equation of motion for Fourier- 
transformed perturbation variable Vk is given by 

c 2 fc 2 _^\ Vk= o. (3. 6 ) 



The above so-called Mukhanov- Sasaki equation (|3.6[) is commonly used to discuss the behavior of the curvature 
perturbation. Here, however, we are interested in the time-dependence of the sound velocity, in which situation the 
potential term z" / z in (|3.6p is not easy to treat since the variable z also contains c s . For example, if we consider a 
step- function-type variation of the sound velocity, a square term of (c' s /c s ) appears in z"/z, which becomes the square 
of the delta function 8(t — To) and makes the analysis impossible where To denotes time when sound velocity changes. 

It is therefore clear that the equation for Vk is not suitable 1 and so we should introduce a new variable Uh, which 
is related to Vk as |47| 

- Cs k 2 u k = z(^)' , Cs v k = e{^)\ (3.7) 



As discussed in Appendix A, this is not the case in the curvaton scenario. 
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where we have defined 

c s z 

The basic equation of motion (|3.6[) for is translated into the new equation in terms of u k , which is 



(3.8) 



u'l + I cik 2 --)u k = 0. (3.9) 

Note that the term (c' s /c s ) does not exist in 6" JQ since the variable 9 does not depend on c s due to the definition of 
9, (|3 . 8|) . We have to solve this equation under the assumption that the background evolution satisfies the slow-roll 
conditions (|2.7p . 

The term 9" /9 is rewritten in terms of slow-roll parameters as 

where we have used slow-roll approximation (|2.7p and a useful equation 

Therefore, we obtain the basic equation as 

u'l + (c 2 s k 2 - u k = 0, (3.12) 



where we have defined 



and approximate it as 



„* = | + e+ I. (3.13) 



IV. STEP-LIKE VARIATION OF THE SOUND SPEED 

In this section, we compute the curvature perturbation for a model such that time variation of the sound velocity 
is described by a step function as 

{ C S 2 (T > T ) 

To take into account the transition of the sound velocity, it is important to impose a matching condition to the solution 
in (|3.9[) at r = tq when the sound velocity suddenly changes. The matching condition is obtained by integrating (|3.9[) 
in an infinitesimaly small time interval [tq — St, tq + St], yielding two conditions; 

Uk(To - St) = Ufe(r + St), u' k (T - St) = u' k (T + St) . (4.2) 

Hereafter, we use the following expression 

Uk(ro - St) u kl , u k (T + St) u k2 . (4.3) 

In the regime when t < tq, setting c s = c s \ leads to the equation of motion 

c 2 sl k 2 - V —^) u k = 0> (4-4) 



and its solution is obtained by 



u k i = \J-kc s iT diH^ (-kc sl r) + d 2 H^ '(-kc sl r) 



(4.5) 



where H^^ 2 \—kc s iT) denote the Hankel functions and rf L2 are constants to be determined by the initial condition 
at r — > — oo. We choose the adiabatic vacuum at the initial time in terms of Vk- 
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From the equation Q3.7p . we find 



well inside the horizon, so we can take 



\/2kc sl 
c sl k 2 u kl = v' k , 



Ukl 



—ikc s \T 



and match it with the limiting form of the Hankel function, 

~2 



x) 



I 2 r / w i\i 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



where the upper and lower signs correspond to the first and the second types, respectively. Hence it leads to the 
choice of d\ , d 2 as 



di 



2k 3 ^ V Csi 



r ( 2v + l ., 

CXp I ; 7TI ) , d 2 = U. 



(4.10) 



Neglecting all the phase factors which is irrelevant for calculating the power spectrum, the solution Uki takes the form 

(4.H) 



= ^^H^(-kc sl r). 



2k 



Next, in the regime when t > tq, setting c s = c s2 leads to 

'a' k H^(-kc s2 T) + P' k HW(-kc s2 T) 



and it is rewritten by 



where we have defined 



u k2 



U k 2 



\/~kc s2 T 



2k 



H^(-kc s2 r) + PkH^(-kc s2 T) 



From (|4.2p . we finally obtain the coefficients a k and j3k as 



Otk 



iixkTQ 
4 

inkro 



c s2 H^{-kc sl T )Hl%(-kc s2 T ) - c sl H^{-kc s2 T )H^l x {-kc sl r ) 

C a 2^ 1) (-*C-l7D)-H'Sl(-*C.27D) - C sl ^ 1 )(-fcc s2 To)^ 1 + ) 1 (-fc Csl T ) 



where we have used the following relation among the Hankel functions, 



(4.12) 
(4.13) 
(4.14) 

(4.15) 
(4.16) 

(4.17) 
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We have to estimate the power spectrum of £ at the final time r — > 0. By using (|3.7j) . the curvature perturbation 
can be rewritten in terms of the variable u as 



z 

a/27 V 2/ 



= u \ u k ~ ~K u k 



(4.18) 



where, at the last transformation, we have chosen the super-horizon limit and neglected u' k term. Therefore, we need 
the final values expressed as u k i(r — s- 0). From the solution (|4.13|) . we obtain 



\/ — 7TT i „, . 



-kc s2 T 



where r(x) denotes the Gamma function and we have used the limiting form of the Hunkcl function as 

H^(y) ~ -H^(y) ~ l -V{u) (|)~" , (y 0). 
By using (|4.18p and (|4.19p . we can calculate the power spectrum of curvature perturbation as 

P c (k) = \C k (r H- 0)| 2 « [aff(l - e)] 2 ^" 1 (fcc*)" 2 " (1 + r?) \a k - p k \ 2 + 0(e 2 , V 2 ), 



(4.19) 



(4.20) 



(4.21) 



where we have used the slow-roll approximation. Hence the term \a k — P k \ 2 is important to determine the power 
spectrum. a k — Pk can be described by using (|4.15[) and (|4.16|) as 



a k - Pk 



(c s 2J v (-kCslTo)J v +l(—kc s 2To) - C s l J v +1 (—kc s lTo) J v (—kc s 2To )J 

+ ^c a2 A r t/(-A:c s iro)J t/+ i(-fcc a 2ro) - c s iAr„ +1 (-fcc s iT ) J„(-A;c S 2To)J 



(4.22) 



where J„, iV^ are the Bessel functions with v = \ + ^ + e. Then we can calculate \a k — /3fc| 2 as 



|a fc -&| 2 = 



tt 2 *; 2 



(c s2 J v (— feCaiTo) Jy + l(-fcc s2 To) - C s i J, y+ i(-fcc s iT )J^(-fcc s 2T )^ 
+ ^c s . 2 A r ^(-fcc s iTo)J 1 , + i(-fcc S 2To) - C sl N v+1 (-kc s lTo)Jv(-kc s 2T )*) 
Here, we define new variables A and ko as 



1 1 

c s2 =Ac s i, -c s2 r = — , -c s1 t = — — , 

fc Afc 



(4.23) 



(4.24) 



and using them leads to 
tt 2 *: 2 



\a k -p k \* = 



\k 2 A 2 



A J y+i T 



■If, 



k 



" \Ak Q J v \Ak 



1^ 



k 

~Aka 



k 
Ah 



v \ko 



k 
Ah 



Ak 
N„ +1 



k 

Ak 



Finally, we can evaluate the power spectrum at the sound horizon crossing aH = A:c S 2 as 



P ( (k) 



2 1 +'?H 



1 t} \ (1 - e)"+ 2e H 2 



r z - + - + e 
1 2 2 



fc 3 



{l + ri)\a k -p k \'' 



(4.25) 



(4.26) 



aH— k c s 2 



In the lowest-order slow- roll approxmation, we can set v = in which case the Bessel functions are expressed by 
trigonometric functions and the modulation factor \a k — p k \ 2 can be recast to the following simple form (see Fig. 1(a)). 



\a k -p k \ 2 = A 



A 2 



- 1 sin z 



(4.27) 
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FIG. 1: (a) The modulation factor in the leading-order slow-roll approximation (|4.27p in the case A = 0.9 (A mp = 0.17). (b) 
CMB temperature anisotropy spectrum for some different values of fco- The other cosmological values are chosen as the WMAP 
7-year mean values. 



In this case, neglecting the slow-roll correction in the numerator, the dimensionless power spectrum (|4.26[) becomes 



Al(k) 



2^ 



H 2 A 

87r 2 ec s s 



1 



1 

A? 



1 sin 2 



(4.28) 



This result recovers the usual power spectrum with constant sound velocity if we take A = 1 . In the large-scale limit 
k — > 0, we can drop the oscillatory term and the final expression becomes 



A 2 (fc) 



H 



8ir 2 ec s i 



which represents the power spectrum for modes which cross the horizon far before the transition time To- 
Taking the result (|4.28p as the input primordial power specrum, or more quantitatively speaking, taking 



-4, 



A 2 (k) = A 

2.43 x 10" 9 , k 



1 4~ Amp sii 
0.002PV1PC- 1 ] 



n s = 0.963, 



(4.29) 



(4.30) 



as the input one, we have computed CMB temperature anisotropy for several values of fco (see Fig. 1(b)). Comparing 
(|4.30p with (|4.28|) . we can clearly see that A mp = 1/A 2 — 1. The parameter fco determines the scale under which the 
rapid modulation appears. As we expect, if the value of fco is chosen as O(0.001[Mpc -1 ]) or smaller, then the CMB 
spectrum starts to oscillate at the relatively large scale around I ~ 10, which may better fit some anomalous data 
points obserbed in WMAP. Actually, as one exapmle value, if we take fco = 0.003[Mpc _1 ], A mp = 0.087 and neglect 
the small scale (£ > 200) oscillation, we have found that the % 2 -value improves 3.2 compared with the usual case of 
power-law primordial power spectrum. 

In Fig. 2, we have plotted the difference between the CMB temperature anisotropy spectrum calculated by the 
primordial power spectrum with oscillations (|4.30[) and the one by the usual power-law power spectrum without 
oscillations, which we denote C? sc and Cf d , respectively. Here the difference is defined as 



AC, 
C, 



c? 



-/</ 



c, 



std 



Note that the overall amplitude A s in C° sc is different from that used in Cf d by a factor (1 + A mp /2)~ 
also plotted the expected scatter of Ce due to the cosmic variance, which is given by 



Ct 



cv 



(2£+l)f sky - 



fsky = 0.65. 



(4.31) 



We have 



(4.32) 
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Note that the observational error of WMAP is cosmic- variance limited up to £ ~ 400 and it will be extended to 
£ ~ 2500 for PLANCK. Hence from this figure we can constrain A mp as A mp < 0.1 now for k = lO-^Mpc" 1 ] which 
corresponds to the largest scale that can be measured with CMB experiments. 




FIG. 2: The difference between CMB temperature anisotropy spectrum from the initial power spectrum with the oscillation and 
without the oscillation. We choose fco = 0.0001 [Mpc ] and the obserbational detection limit in the PLANCK-like experiment 
is also depicted. 



This feature in the spectrum can be compared with the Trans-Planckian signatures [22h24| . with the particle 
production effect, or with the spectrum generated in the case the inflaton potential has a break in its slope. In the 
first case, the modification of the initial vacuum can change the relation and then leads to the ringing signatures in 
the primordial power spectrum. The exact form of the pattern depends on the choice of the hypersurface on which 
the initial condition is imposed. In the so-called Boundary Effective Field Theory (BEFT) approach |48|,|49|, the final 
expression becomes 



7W T (fc) = A s 



sin 2- 



a.iA V QiH,, 



(4.33) 



where a, and Hi are the scale factor and the Hubble parameter on the initial condition hypersurface, respectively, and 
A is a cutoff scale. This is similar to our result (|4.28p , but the amplitude of the oscillation depends on the wavenumber 
k, which makes crucial difference with our result. 

In the second case, Barnaby et al. [26[ discuss that the particle production (PP) make a bumplike contribution to 
the power spectrum, which is well fitted by 

k \ ™ s_1 / k \ 3 / irk 2 



rMk)=A s [-) + A^-J exp^j. (4.34) 

The feature is localized around fc*, which makes sharp difference with our result. 

In the last category, where the second derivative of the inflaton field has a sudden change at t = tq , the modulation 
factor can be calculated similarly and the final expression is also similar to our result at the first glance. However, 
there is a big difference because in this case the slow-roll parameter has a step-like function over the transition time 
and the tilt of the spectrum takes the different values before and after the transition [H, H3] . 
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TOP-HAT TYPE VARIATION 



In this section, we consider the other type of variation of sound velocity, namely, the top-hat type. This is a simple 
extension of the model discussed in the last section. We parametrize the variation as 



Cgl (T < T ) 

C S 2 (To < T < To) . 
C S 1 K < T ) 



(5.1) 



and impose the same matching condition to the wavefunction Uk as in (|4.2j) at the two transition times To and Tq. 
The solutions of the equation of motion during constant sound velocity regimes are expressed in the same way as in 
the last section and we parametrize them as follows; 



P^Q^HM(-kc sl T) (r<r ), 



2A- 



Uk2 



Uk3 



2k 
2k 



a kl H^\-kc s2 r) + (i k iH^\-kc s2 T)\ (r < r < r ), 
a k2 H^(-kc sl T) + P k 2H^(-kc sl T)] (r < r), 



(5.2) 
(5.3) 
(5.4) 



where we have chosen the Bunch-Davies vacuum state at r <C tq. Two coefficients a k i and (3ki hi (|5.3p correspond to 
a k and j3 k in (|4.13l) . Since the solutions and the matching conditions are the same, the coefficients a k \ and f3 k i are 
given by (|4.15|) and (|4.16j) . The matching conditions at r = Tq gives the relations between a k 2,Pk2 and a k i, f3 k i', 



Otk2 



c s2 H^l 1 (-kc S 2T l )H { J- ) {-kc s iT' ) - c s iH { ^l 1 (-kc sl T' )Hl 1 \-kc s2 T' ) \a k i 



+ CaH^i-kc^H^i-kCsni) - c sl Hl+i(-kc s i4)Hl 2 \-kc s2 4) 



(5.5) 



inki 



1,2 



c s2 ^+ ) 1 (-fcc s2 r )i7( 1 )(-fcc sl r ) - c sl Hi%(-kc sl T^H^(-kc s2 T^ \a kl 



(5.6) 



We want to know the power spectrum of the curvature perturbation at the final epoch r — >■ 0. As we have already 
seen, all we need is Uk3 (r — ^ 0). Combination of the solution (|5.4[) and the limiting form of the Hankcl function leads 
to 



u k 3{T -> o) ~ — — r(i/) 



2k TT ' ' V 

and the power spectrum of the curvature perturbation is described as 



(a k2 - Pk2), 



V\ 2 



P c (fc)^|C(r^0)| 2 = |-(l + ^) 



(5.7) 



(5.8) 



Also in this case, the factor to determine the modulation pattern is \a k2 — Pk2\ 2 , which can be calculated by inserting 
(|4.15l) and (|4.16p into (|5.5[) and (|5.6[) . The final functional form however is so complicated. We therefore show the 
result only in the leading order slow-roll approximaton; 

1 



A 2 sin 2 (-fcc s iTQ) + cos 2 (-fcc s iTo; 



\a k2 - p k2 \ 2 ^ cos 2 [-kc s2 (ro - r )] + sin 2 [-fec s2 (ro - r ) 

-2^A- -^J cos(-fcc s iro)sin(-fcc s iro)cos [-fcc s2 (ro - r )] sin[-fcc s2 (ro - r )] + £>(e,r?). (5.9) 



As expected, the modulation is determined by the combination of sine and cosine osillations. We can transform 
5.9[) to the following form. 
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FIG. 3: (a) The modulation factor in the leading order slow-roll approxmation (|5.10|l in the case A = 1.1 and B = 0.9. (b) 
CMB temperature anisotropy spectrum for some different values of fco- The other cosmological values are chosen as the WMAP 
7-year mean values. 



a 2 - 1 r 

\a k2 - Pk2? =1 - -^j- (A 2 - l)sin 2 [-kc s2 (^ - r )] - {A 2 + 1) sin 2 (-fc Cs iT ) 

(A + l) 2 (A — l) 2 
H x sin 2 [-kc s ir' Q + fec s2 (ro - r )] H sin 2 [-fcc gl To - kc^r^ - r )] 



(5.10) 



It is clear in this expression that there is no modulation when A = 1 or To = Tq, which means no change in 
sound velocity. If we consider the long- wavelength limit (fc — > 0), the dimensionlcss power spectrum of the curvature 
perturbation reduces to the same form as (|4.29p in the lowest-order slow-roll approximation of the numerator, which 
is valid because the long-wavelength mode exits the Hubble horizon when the sound velocity is c s \ and superhorizen 
mode does not feel the sound velocity change. 

The modulation factor (|5.9p or (|5.10|) is plotted in Fig. 3(a) and the CMB temperature anisotropy spectrum is 
plotted in Fig. 3(b). In both figures, we defined three parameters ko, A and B as — c s ir' a = 1/ko, To = Btq and 
c S 2 = Ac s i, respectively. In these figures, We can see a beat, which is very natural because we now have two 
charactarizing parameters c s \Tq and c S 2(tq — To) that determines the ocsillation periods and the final modulation 
pattern is the superposition of these waves. This result is easily extended to more general cases; if the variation of the 
sound velocity occurs suddenly and during it remains constant at other regimes, then the modulation of the power 
spectrum is some superpositions of several waves composed by trigonometric functions, which appears as a beat and 
the amplitude of the oscillation is independent of the wavenumber. 



VI. SUMMARY AND DISCUSSION 



In summary, we have calculated the power spectrum of the curvature perturbation in case the sound velocity 
changes suddenly during inflation. In order to focus on the effects of the change of the sound velocity, we have 
assumed that the slow-roll parameters that reflect the background evolution are not affected and remains constant 
during the transition of the sound velocity and found that in the shorter-wavelength modes (with larger k) compared 
with the Hubble radius at the transition time there appears osillation patterns due to the mode mixing. In the 
simplest model in which the sound velocity experiences only one transition, the oscillation in the power spectrum is 

expressed as oc 1 + A mp sin 2 where k corresponds to the wavenumber that crosses the sound horizon at the 

transition time, A mp is defined as A mp ~ I /A 2 — 1 and A is the ratio of sound velocity before and after the transition. 
Different from the trans-Planckian signatures or the models with a sudden change of the slow-roll parameters, the 
amplitude of the oscillation does not depend on the wavenumber, which gives a clue to find an evidence of the change 
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A— n 995 
(A mv = 0.0099) 


A=0.99 
(A m „ = 0.0197) 


A=0.97 
(A™, = 0.0574) 


fco = 0.0001 


-0.054 


-0.11 


-0.27 


k = 0.0003 


-0.063 


-0.11 


-0.13 


k = 0.001 


-0.34 


0.60 


-0.81 


fco = 0.003 


0.045 


0.30 


3.34 



TABLE I: This table shows the Ax 2 values in the step- function model described in Sec lIVI compared with WMAP 7- year best 
fit cosmological parameters in the simple ACDM model. The unit of the wavenumber fco is Mpc -1 . For the calculation of the 
X 2 values, we exploited the WMAP-7year likelihood funcion [5C( |. 





A = 1.005 


A = 1.01 


A= 1.03 


k' = 0.00001 


-2.26 


-1.54 


2.53 


k'o = 0.00005 


-1.88 


-1.40 


1.28 


k'o = 0.0001 


-1.72 


-1.25 


1.33 


k'o = 0.0005 


-2.01 


-1.59 


2.28 



TABLE II: This table shows the A\ 2 values in the top-hat function model described in Sec|V] compared with WMAP 7- year 
best fit cosmological parameters in the simple ACDM model. We have defined a new variable k' as — c s \To = l/k' . The 
unit of the wavenumber k' is Mpc -1 . Here, we fix fco = 0.001[Mpc -1 ]. For the calculation of the \ 2 values, we exploited the 
WMAP-7year likelihood funcion [H(|. 

of the sound velocity by some observations. Actually, we have computed the CMB temperature anisotropy using the 
primordial power spectrum we have obtained and we can see the oscillating patterns in the resulting angular power 
spectrum originating in the modulation of the primordial curvature perturbations. 

These oscillatory behaviors make it difficult to constrain the amplitude of the modulation using Markov Chain 
Monte Carlo (MCMC) analysis due to the problem of the convergence. Thus, we calculated the difference between 
the CMB temperature spectrum from the oscillatory primordial power spectrum and the one without the oscillation 
and compared it with the scatter due to the cosmic variance. We have confirmed that if ko = 10 -4 [Mpc -1 ] the current 
CMB experiment is sensitive to A mp > 0.1 or in other words is sensitive to A < 0.95. 

To discuss how WMAP experiment can constrain or find the amplitude of the oscillation of the primordial power 
spectrum, in Tables I and II, we have computed A\ 2 for various parameter values in the two models we have 
considered. We can check from this tabic that in the step-function model, too small a transition scale (or too large a 
transition wavenumber fco) and too large an oscillation amplitude cannot fit the present CMB anisotropy data. In the 
top-hat function model, we have found several good improvements of Ax 2 values, though we have more parameters 
than step-function model in thic model. Although we have been unable to find parameter values that satisfy Akaike's 
Information Criterion (Hlj . several negative Ax 2 values may give us some motivations to search for the parameter 
values for much better fit to WMAP data or more accurate future CMB experiment data in our model. 
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Appendix A 

In this appendix, we comment on the consistent background evolution which keeps slow-roll condition during the 
sudden transition of the sound velocity. In general, the slow-roll conditions are violated by changing the sound velocity 
so rapidly like step-type funtion without numerical fine-tuning between non-canonical kinetic terms. 

One way to overcome the difficulty is to consider the curvaton scenario [52l. [53j . If the curvaton Lagrangian has the 
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following form, 

C = \d^d^ + \f(4>)d^a + £7i(0)(0„ W) 2 + \g2{4>){d^o) 2 

1 1 771 2 

+ -g 3 (0)(d^d^)(d„ad»a) + -g 4 (<P)(d^a) 2 + ■■■- V(4>) - —a 2 , (6.1) 

where <f> and a denote the inflaton and the curvaton, respectively, and gi{4>) = 0(A~ 4 ). Assuming that, during 
inflation, the curvaton is set to some constant value, we can write down the second-order Lagrangian of the curvaton 
perturbation 

£ ( 2) = £W! ff , 2 _ JW + GW 2 _ ^ 2 

2a 2 2a 2 2 

Here, 

F (<f>) 2 = f(<f>) + (gM + .94(0)) X+ + ---, (6.3) 
G(0)=-54(0)^ + ---, (6.4) 

where = <j> 2 /2. Then we can read off the sound velocity for the curvaton perturbation as 

2 f(0) 2 + c(0) m+gmx* (c - 

F(^) 2 m + faM+gAiMXi [ - } 

In order for c s to deviate from unity, there should be difference between the coefficients of time-derivative and of 
space-derivative. During inflation, the background breaks the symmetry between time and space through <fi. Because 
of the coupling d^tpd^a, the sound velocity of a can be affected by (f>. Actually the important factor to cause the 
deviation from unity in (|6.5[) is 54(0), namely, the coefficient of the coupling d^ipd^a. This observation also implies 
that it is difficult to change the sound velocity without affecting background evolution in generic single field inflation 
models since, in such models, all the higher order couplings contribute to both time and space derivatives of the 
inflaton perturbation and hence the deviation in the sound speed is related to a form of a background equation of 
motion. 

To be consistent with the assumption that the background value of the curvaton keeps constant, the condition 

o 1 der 

IT* = ~,dN « 1 (6 - 6) 

should be satisfied, where N stands for the number of e-folds and an over dot denotes differentiation with respect to 
the physical time. The background equation of motion for a is 

d 2 er ,„ „ x da / m \ 2 

+ (3-e + 2a)— + ( — a = 0, (6.7) 



dN 2 y 'dN \HFJ 



where a = F/HF. Imposing <C 1, this equation tells 

da 1 



dN 3-e + 2a \HFJ 

Thus the condition (|6.6[) is satisfied if m/H <C F. 

As for the curvaton perturbation, the second order Lagrangian can be written as 



f^) 2 - (6-8) 



~gC^ = a 2 F 2 



l „' 2 1 Jtv„\* m a - 2 
2° -2° s(W) 2F 2 " 



(6.9) 



Introducing the new variable v a = z a a where z a = aF which corresponds to v for the inflaton scenario in the main 
text, (|6.9[) can be rewritten as 
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Now we can derive the evolution equation for the curvaton perturbation in Fourier space as 



+ I c 2 k 2 + 



2F 2 



Vak = 0. 



(6.11) 



Imposing the above-derived condition m/H -C F, we can neglect the mass term in this equation during slow-roll 
inflation. On the other hand, the second derivative of z a can be evalutated as 



f = ( aH Y 

Zfj 



da 

(1 + a)(2- e + a) + — 



5.12) 



If we want to change the value of sound velocity without affecting the background evolution, F(<p) should be fixed. 
This is achieved when /(</>) = 1 and 33(0) = —gi(4>) = g ((/>)■ In this case, c 2 = 1 + g^X^ + ■ ■ ■ and we can 
change only c s by choosing g(cj>) as an appropriate form with F = 1. Actually, now that z'l/z a = {aH) 2 (2 - e), the 
techniques to calculate the power spectrum introduced in the main text can be easily applied without introducing the 
new perturbation variable Uk which is no longer convinient to use. 
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